Numerical studies of the PfafRan model of the = | fractional quantum Hall effect 
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The PfafBan model has been proposed for the fractional quantum Hall effect (FQHE) at — 
|. We examine it for the quasihole excitations by comparison with exact diagonalization results. 
Specifically, we consider the structure of the low-energy spectrum, accuracy of the microscopic wave 
functions, particle-hole symmetry, splitting of the degeneracies, and off-diagonal long range order. 
We also review how the 5/2 FQHE can be understood without appealing to the PfafRan model. 
Implications for nonabelian braiding statistics will be mentioned. 



I. INTRODUCTION 

5/2 is the only even denominator fraction securely ob- 
served in a single layer fractional quantum Hall effect 



(FQHE 



)b2,^ 



(The fraction 7/2 is trivially related to it by 



particle hole symmetry in the second Landau level.) The 
model of noninteracting composite fermions predicts a 
compressible Fermi sea at half filled lowest Landau level, 
which provides a good description of the compressible 
state here^. A promising scenario for the incompress- 
ible state at the half filled second Landau level at = 5/2 
is based on the idea of pairing of composite fermions, de- 
scribed by a Pfaffian wave function^ii. Several studies 
have supported this interpretatio n^i^d^d^d^ . 

The Pfaffian wave function is the exact ground state 
of a singular three-body model interaction (cf. Eq. [2] be- 
low). Exact solutions for quasiholes are also available for 
this model interaction. A case has been made, both from 
analytical argument o^d'^d^ and numerical calculations^^, 
that these Pfaffian quasiholes have the remarkable prop- 
erty of nonabelian braiding statistics. Recently, the non- 
abelian statistics has taken additional importance be- 
cause of proposals to test it experimentall y ^^d^d^ , and 
to exploit it for quantum computatio n^^d^i^Pi^^ . That 
makes it important to perform an examination of the 
applicability of the Pfaffian model to the real. Coulomb 
solution. The Coulomb ground state wave function has 
been compared to the Pfaffian ground state wave function 
in the past-'^^ and found to have overlaps in the range 
0.69-0.87 for 8 to 16 particles. Our recent comparisons of 
the Pfaffian quasiholes and the real Coulomb quasiholes^^i 
showed a worse agreement. Because the route to non- 
Abelions is via the Pfaffian model and the degeneracies it 
implies, these studies have relevance to nonabelian braid- 
ing statistics as well. 

This paper briefly reviews our previous work, at the 
same time providing many new results relevant to this 
problem. In Sec. [TTl the Pfaffian model is defined and 
some relevant results are reviewed. In Sec. IIIII we com- 
ment on the particle-hole symmetry violation by the Pfaf- 
fian family of states. In Sec. IIVI the absence of off- 
diagonal long range order in the Pfaffian state is pointed 
out, and the relevance of this finding is discussed. In 
Sec. I VII we check the assertion, commonly made in the 
literature, that the energy difference of the quasihole 



states, that are degenerate for the three-body model in- 
teraction, remains exponentially small for the Coulomb 
interaction. In Sec. IVIII we attempt to separate the pos- 
tulated charge i quasiholes for the model interaction as 
well as the Coulomb interaction. Finally, in Sec. IIXI we 
elaborate an alternative approach for the explanation of 
the I FQHE. Short reports on parts of this paper have 
appeared elsewhere^i^^. 



II. THE PFAFFIAN MODEL 

Throughout this article we will assume that the low- 
est Landau level (LL) is full and inert, and the two- 
dimensional electron gas in the second LL is fully po- 
larized. All calculations are performed in the spherical 
geometry. The objective is to determine the ground state 
and the low-energy excitations for the Coulomb interac- 
tion 
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in second LL at filling factor j/ = 1/2 (e is the static di- 
electric constant of the host semiconductor). This prob- 
lem is equivalent to electrons in the lowest Landau level 
interacting with an effective interaction . We will use 
the lowest LL to simulate the second LL physics in what 
follows. 

The "Pfaffian model" considers a three-body model 
interaction^'14^ which in the spherical geometry takes the 
form 
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where PijkiLmax) is the projection operator onto an elec- 
tron triplet with orbital angular momentum Lmax = 
3Q — 3. The angular momentum Lmax corresponds to 
the closest possible configuration of an electron triplet. 
Thus, V^^ does not penalize the closest approach of two 
electrons, but there is an energy cost when three electrons 
are in their closest configuration. 

This model has a unique, zero energy ground state at 
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1/2 (Moore and Read*'): 



III. PARTICLE-HOLE SYMMETRY 



Pf 



<&! = l[{z^ - Z,), (3) 



where "Pf refers to "PfafRan." This wave function de- 
scribes a paired state of composite fermions. This model 
also produces exact zero energy eigenfunctions for quasi- 
holes, as the flux through the sphere is increased. These 
zero energy eigenstates are referred to as the "PfafRan 
quasihole (PfQH) states." As [V^^,L^] = 0, the states 
spanning the PfQH sector may be chosen with a definite 
orbital angular momentum L. (See Ref. for a thorough 
study of the PfQH sector on the sphere.) Appropriate 
linear combinations of these states produce spatially lo- 
calized quasiholes. For two quasiholes at rji and 772, the 
wave function is given by^ 



Pf 

2-qh 



= Pf 



{zj ~ r]i){zj - ri2) + (z ^ i) 

(Z, - Zj) 



(4) 



For two coincident quasiholes, rji ~ ri2 = Vi this reduces 
to a charge ^ vortex: 



(5) 



Separately, each quasihole has a charge deficiency of j 
associated with it. Unlike for the vortex, the density 
does not vanish at the center of a quasihole. Analogous 
wave functions can be written for an even number (2m) 
of quasiholes. Exact wave functions for quasiparticles are 
not available. 

Several wave functions can be associated for a given 
configuration of 2m quasiholes, which correspond, in the 
appropriate generalization of Eq. (|4]) to 2m quasiholes, to 
different ways of grouping half of the rj^. 's with z^ and the 
other half with Zj. It has been shown^^ that only 2™"^ 
of these functions are linearly independent. Adiabatic 
braiding of quasiholes (which is feasible for a gapped sys- 
tem) can take the system from one linear combinations 
of PfQH states to another, which lies at the origin of 
nonabelian statistics of quasiholes. 

To study bulk properties, it is convenient to formu- 
late the problem of interacting electrons in the spheri- 
cal geometry, in which the electrons move on the surface 
of a sphere and a radial magnetic field is produced by 
a magnetic monopole of strength Q at the center.— 
Here 2Q<j)o is the magnetic flux through the surface of the 
sphere; (po — hc/e, and 2Q is an integer by Dirac's quan- 
tization condition. Then wave functions in Eqs. (I3][5]), 
which are written for the disk geometry, can be mapped 
to the sphere by the stereographic mapping^^, which 
amounts to the substitution 

{Za - Zb) ~* {UaVb ~ VaUb), (6) 

for all coordinate differences, where Ua — cos ^e"*"^"/^ 
and Va — sin^e*'^''/^ ^j-e spinorial coordinates on the 
sphere. The orbital angular momentum quantum number 
is denoted by L. 



The exact Coulomb eigenstates in any given Landau 
level satisfy particle-hole symmetry, i.e., the exact eigen- 
states at v and l — v are related by particle- hole transfor- 
mation. The wave functions in the CF theory^^ satisfy 
particle hole symmetry to a very good approximation, 
even though there is no symmetry principle that so re- 
quires. For example, the wave functions &i v — n/ {2n — 
1), given by ^'„/(2„-i) = 'PLLL^'?[^'n]*, are almost identi- 
cal to the those obtained by particle-hole transformation 
of the wave functions 5'„'/(2n'-i-i) — '^LLL'&i^'n', with 
n = n' + 1. 

The three body interaction does not satisfy particle- 
hole (p-h) symmetry. To get a feel for the extent to 
which this symmetry is broken, we have considered the 
system of TV = 8 particles at 2Q = 15. In this case, 
particle hole transformation gives eight holes (to be dis- 
tinguished from quasiholes) at 2Q = 15. We obtain the 
exact spectrum of the V^^ model interaction, which is 
given in the upper left panel of Fig. [TJ This system cor- 
responds to four quasiholes, and has a number of zero 
energy states, which form the Pfaffian quasihole sector. 
We obtain the particle-hole conjugate of each eigenstate, 
called and calculate its energy expectation value for 
the V^^ interaction. When there are several degenerate 
Pfaffian quasihole states, we diagonalize V^^ in the sub- 
space of the p-h conjugate states to obtain the energies. 
The resulting spectrum is shown in the top right column 
of Fig. [TJ (For the Coulomb interaction, this exercise 
would produce a spectrum identical to the original one, 
apart from an overall energy shift.) We construct sym- 
metrized states cx (^ -|- 'I'^), which satisfy particle- 
hole symmetry by construction; the resulting spectrum 
for these states is given in the lower left panel of Fig. [T] 
The lower right spectrum is for antisymmetrized states 
^fa oc (* - '^''). 

Tabic [1] shows the squared overlaps between the orig- 
inal Pfaffian quasihole states with the various states ob- 
tained with the help of p-h conjugation. To handle the 
multiplicity of the PfQH sector for L = 0,2,4,6 (cf. 
Fig. [T]), the overlap between two subspaces has been de- 
fined in a basis-independent manner (see caption of Table 
[ij. The overlaps are not particularly high; for example, 
in the L = Q part of the quasihole branch, which contains 
two states for = 8, the overlap is 0.511, and deterio- 
rates for higher L's. Similar numbers are obtained for 
other states in the PfQH sector. The near orthogonality 
of 5* and ^E*^ at L = 8 is accompanied by a very high 
energy of ^E''^. 

These results demonstrate a substantial breakdown of 
the p-h symmetry by the V^^ interaction. The Pfaf- 
fian quasihole band is absent in all of the new spectra; 
the states derived from the Pfaffian quasihole band are 
mixed up with other states. It has been shown^i that the 
particle-hole symmetrization {'^ — » '^^) of the Pfaffian 
wave function improves the overlap with the Coulomb 
ground state. Our results show, however, that this also 
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destroys the degeneracy of PfQH sector. One can ask 
whether the nonabehan statistics of the Pfaffian quasi- 
holes is robust to p-h symmetrization; we are not able 
resolve this question definitively by a direct calculation 
of the braiding phases, which requires much larger sys- 
tems. 
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posite fermions opens a gap to produce FQHE but does 
not establish long range phase coherence in the electronic 
state. For this purpose, we calculate the off-diagonal long 
range order parameter: 

where ?/'(r) and ip^{r) are the usual annihilation and cre- 
ation field operators. We place the primed coordinates 
near the north pole, separated by a distance equal to 
the magnetic length, and the unprimed coordinates at 
the south pole, also separated by a distance equal to the 
magnetic length. The results in Table |TT1 obtained by 
Monte Carlo calculation, demonstrate the absence of off- 
diagonal long-range order in the Pfaffian wave function. 
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FIG. 1: Upper left panel shows the original spectrum of V 
with iV = 8, 2Q = 2iV - 1 (four quasiholes); the Pfaffian 
quasihole states have zero energy. Also shown are the spectra 
for the p-h conjugate states (top right), the p-h symmetrized 
states (bottom left), and p-h antisymmetrized states (bottom 
right). The diamonds show the states derived from the Pfaf- 
fian quasihole branch. 
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TABLE I: Squared overlaps between the subspaces spanned 
by the zero-energy states and the subspaces spanned by 
their particle-hole conjugate, particle-hole symmetrized, and 
particle-hole antisymmetrized images, respectively. Squared 
overlaps are defined as O = K*4-qh,i|*4_qh,j)P/A/', 
where A/" is the number of degenerate multipletsi^ of V^^ at 
L, and i, j — 1, ■ ■ ■ ,J\f. 



IV. OFF-DIAGONAL LONG RANGE ORDER 

We wish to stress that the Pfaffian wave function does 
not represent a true superconductor; the pairing of com- 



TABLE II: Off-diagonal long-range order parameter 
G(ri, r2,ri, r2) with ri and r2 separated by Ib about 
the north pole, and r'l and r'2 separated by Ib about the 
south pole for the paired CF wave function Pf(l/(2:i — Zj))$i. 



V. TESTING THE PFAFFIAN QUASIHOLE 
WAVE FUNCTION 

We have recently carried out comparisons between 
the Pfaffian and Coulomb quasiholes^^. Figs. [21 [3l 
and [4] show the spectra for states with two and four 
quasiholes for TV = 10 and 12 electrons. For 10 elec- 
trons, the Pfafhan model predicts zero energy states at 
L = 1,3,5 and L = 0^, 1°, 2^, 3\ 4^, 5^, 6^, 7^, 8^, 9°, 10\ 
respectively (the superscript denotes the degeneracy) , for 
two and four quasiholes. These states form the Pfaffian 
quasihole band. For 12 electrons, the Pfaffian quasihole 
band contains states at L = 0,2,4,6 for two quasiholes 
and L = 0^, fO, 2^, 32,45, 5^, 6^ 7^, 8^,91, 10^, 11°, 12^ for 
four quasiholes. For 14 electrons, the Pfafhan quasihole 
band for two quasiholes has states at angular momenta 
L = 1,3,5,7. 

The Coulomb spectra in Figs. [2J3] do not show 
well defined bands that have a one-to-one correspon- 
dence with the Pfafhan quasihole bands. Ref. [22 gives 
overlaps between the Pfaffian and Coulomb quasihole 
states, which are generally worse than for the ground 
state. Fig. [5] depicts for the two quasihole state 
(for 14 electrons) the "total overlap," defined as O = 
EL=i,3,5.7l(*2-qhl*™ui>lV4 where Ivf^qh) is the two 
quasihole state with — L and l^'^oui) ^^1*^ Coulomb 
ground state with — L. This figure shows the depen- 
dence of the overlap on the form of the interaction; by 
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FIG. 2: Spectra at = | for the model interaction (left 
column) , and the Coulomb interaction (right column) for TV = 
10 particles at 21 = 18 (top row) and 21 = 19 (bottom row). 
For the V^^ interaction, two (four) quasiholes are expected 
for 21 = 18 {21 — 19). The spectra on the left were also given 
in Ref. 12. This figure is taken from Ref. HI. 
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FIG. 3: Spectra at = | for the model interaction V^^ (left 
column), and the second Landau level Coulomb interaction 
(right column) for A'^ = 12 particles for two (upper row) and 
four (lower row) quasiholes. 



increasing the Vi pseudopotential of the Coulomb inter- 
action by 0.03 units, it is possible to increase the overlap 
from 0.3 to 0.6. For large 6Vi, the solution is essen- 
tially the lowest-LL Coulomb solution; Fig. [S]thus shows 
that the Pfaffian wave functions provide a comparable 
description of the state in the lowest two Landau levels. 
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FIG. 4: Spectra a.t u — ^ for the model interaction V (left 
column), and the second Landau level Coulomb interaction 
(right column) for A'' = 14 particles for two quasiholes. 




FIG. 5: The overlap between the low-energy excitations of the 
second-LL Coulomb and V^' interactions for A*' = 14 particles 
as the leading pseudopotential Vi is changed. 



VI. ENERGY SPLITTING OF THE PFAFFIAN 
QUASIHOLE STATES 

The Pfaffian model predicts a 2™"-'^ degenerate wave 
functions for any given configuration of 2m quasiholes, 
which is responsible for the emergence of nonabelian 
braiding statistics. Any deviation from the model inter- 
action V^^ lifts this degeneracy, but a case can be made 
that if the energy splitting of these states remains expo- 
nentially small as a function of the distance between the 
quasiholes, the idea of nonabelian statistics remains ex- 
perimentally relevant. It would be of interest to test how 
the splitting behaves in a realistic calculation. Unfortu- 
nately, a good model for the Coulomb quasiholes is not 
available, and it is not known how separated quasiholes 
can be produced in exact diagonalization studies- ^ (Sec- 
tion VII). We study how the Coulomb interaction splits 
the degeneracy while restricting to the PfQH sector. In 
light of the above comparisons, such a restriction is not 
necessarily a valid approximation, because the Coulomb 
interaction causes a substantial mixing with states out- 
side of the PfQH sector. However, a more accurate cal- 
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culation is currently not feasible. 

The calculation requires at least four quasiholes, which 
we place on the sphere at maximal separation, i.e. at the 
vertices of a regular tetrahedron. The Coulomb interac- 
tion in the first and second LLs is diagonalized in the 
space spanned by two Pfaffian quasihole wave functions. 
The overlap and interaction matrices are calculated by 
Monte Carlo methods; an orthonormal basis is found by 
the standard Gram-Schmidt procedure; and the interac- 
tion is diagonalized in this basis. The Coulomb interac- 
tion in the second LL is simulated in the lowest LL by 
an effective interaction of the form 
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where the coefficients Ci are fixed so that the lowest LL 
pseudopotentials^'^ of V°^{r) reproduce all of the second 
LL Coulomb pseudopotentials vSP for odd integral val- 
ues of TO. (For relevant formulas, see Ref. H^.) 

As apparent in Fig. [51 the lowest LL Coulomb interac- 
tion and the effective second LL interaction give different 
results for small {N < 30) systems. Because the en- 
ergy splittings are very close in the 30 < < 54 range, 
we study larger systems {N > 54) with the lowest LL 
Coulomb interaction only. It is likely that the long dis- 
tance behavior of the splitting does not depend on the 
Landau level index (given that the interaction at long 
distances is independent of the LL index) . Fig. [7] shows 
the lowest LL splitting. 
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FIG. 6: The energy splitting of the two four-quasihole wave 
functions for Coulomb interaction. 

The energy splitting is a nonmonotonic function of N 
(or R). Near the local minima the error in the logarithm 
of the energy splitting is seen to become very large. We 
therefore ask how the value of the ODLRO parameter 
at the local maxima decays with distance. While incon- 
clusive, our results are most consistent with a power law 
decay of the splitting: A straight line fits at all the four 



FIG. 7: Same as Fig.|6l with both scales logarithmic. See line 
fitting on Figure [8l 
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FIG. 8: Line fitting on the log-log graph of the energy splitting 
as a function of the distance. A straight line fitted on the local 
maxima of the data is consistent with a power law decay with 
exponent a = —2.37(6). 



bumps in the log-log plot (Fig. [5]) , but not in the semilog 
plot (not shown). A study of larger numbers of parti- 
cles will be required for further confirmation, which is 
impractical at this stage, but assuming a power law, the 
energy splitting decays with an exponent a — —2.37(6). 
We stress again that the fact that the Coulomb inter- 
action causes a substantial mixing with the non-PfQH 
sector diminishes the value of the calculation presented 
in this section. 



VII. SEPARATING QUASIHOLES 

For the purpose of braiding statistics it is necessary to 
consider spatially localized states of quasiholes. In Ref. 
[22I we have studied states of two quasiholes in the pres- 
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ence of delta function impurities that attract the quasi- 
holes. We take the impurities to be placed at one or both 
of the poles so the eigenfunctions have a well defined 
(although they do not have a well defined L quantum 
number). We also assume sufficiently weak strengths for 
the impurity potential, so they do not cause a mixing of 
the Pfafhan quasihole states with higher energy states. 
Our principal results are as follows. 




o.s ' ' ' ' ' ' ' — ' ' — ' — ' — ' — ' — ' — ^ 

0.5 1 1.5 2 2.5 3 0.5 1 1.5 2 2.5 3 



e (rad) e (rad) 

FIG. 9: Left panel: Charge densities of two quasihole states 
for (a) the PfafEan wave function with two quasiholes at two 
poles; (b) the ground state of V^^' with two delta function im- 
purities; (c) the Coulomb ground state with one delta func- 
tion impurity; (d) the Coulomb ground state with two im- 
purities. The impurities are placed on the two poles (or one 
pole in case of a single impurity), so the eigenstates have a 
well defined Lz- The results are for A'' = 12 (dashed lines), 
= 14 (dotted lines) and A*" = 16 (solid lines) electrons at 
2Q — 2N — 2. The density in (a) is calculated by Monte 
Carlo, and in other panels from exact diagonalizaion. When 
the ground state has Lz ^ 0, there are two degenerate states 
at ±Lz; we have shown only one of them for simplicity. The 
normalization is chosen to ensure that the integrated density 
equals A^. Right panel: the integrated excess charge for each 
density, normalized so that the total charge excess is i. Two 
spatially separated quasiholes will exhibit a step at charge \ , 
as approximately seen in the top panel. This figure is taken 
from Ref. IH. 

For the V^^ model a single delta function impurity in 
the lowest LL localizes a vortex (which is a combination 
of two quasiholes) rather than a single Pfaffian quasihole 
for the following reason. The energy of a given wave 
function is equal to a properly weighted average of the 



densities at the positions of the delta functions (for weak 
impurity strengths). For a delta function at [U,V), the 
lowest energy state (which has zero energy independent 
of the strength of the delta impurity) is the one in which 
both quasiholes localize at (J7, V), producing a vortex 
with vanishing density at {U, V). Surprisingly, as seen in 
Fig. m^b), even two delta impurities fail to separate two 
quasiholes, even though the systems are sufficiently large 
at least for the charge- 1/4 Pfaffian quasiholes to be well 
separated (top panels) 

For two Coulomb quasiholes, at first sight, one may 
expect that even a single delta function should produce 
well separated quasiholes, because it can bind one of 
them, which then should repel the other. As seen in 
Fig- [He, d), neither one nor two delta functions produce 
separated quasiholes. In fact, the charge profile is prac- 
tically identical for the two cases. The situation is more 
restrictive for the Coulomb interaction because, instead 
of many degenerate states, we have a single ground state 
multiplet with a definite L. All that weak disorder can do 
is cause a mixing between the different Lz components 
of the ground state multiplet. For the case of two delta 
functions at the two poles, L^ is a good quantum num- 
ber, so the delta functions only lift the degeneracy of the 
Lz states. The lack of quasihole separation in space is 
attributable to the fact that the ground state now has a 
more or less definite L. The absence of exact degeneracy 
inhibits quasihole localization. 



VIII. IMPLICATIONS FOR BRAIDING 
STATISTICS 



The PfafRan quasiholes are believed to obey non- 
abelian braiding statistics. Our finite system studies of 
the Coulomb solutions do not provide a clear confirma- 
tion of the PfafRan model, and therefore of the non- 
abelian statistics. We cannot rule out the possibility 
that the PfafRan physics will be recovered in the ther- 
modynamic limit. It is useful to recall, in this context, 
how the fractional abelian braiding statistic s2Z . ij^. i2P^ of the 
quasiparticles of the v = nj {2n + 1) states has been con- 
firmed theoretically. There, the CF theory provides a 
qualitatively valid description the quasiparticle band, as 
well as accurate wave functions. These wave functions 
are then used for large systems to establish the abelian 
statistics^-'^^. These calculations also demonstrate that 
the braiding statistics is not well defined when quasi- 
particles are overlapping, which is why its confirmation 
requires large systems. The nonavailability of accurate 
wave functions for the 5/2 quasiparticles or quasiholes 
prevents similar calculations of their braiding properties. 
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IX. AN ALTERNATIVE APPROACH FOR 5/2 
FQHE 

It is not known how the Pfaffian wave functions can 
be improved for the two body Coulomb interaction, due 
to lack of variational parameters. Further, the pairing 
of composite fermions is viewed as arising from an in- 
stability of the CF Fermi sea^iiSi^, but the CF Fermi 
sea is not a limiting case of the Pfaffian wave function. 
These observations have motivated us to approach the 
5/2 FQHE from the CF Fermi sea end, without assuming 
any pairing at the outset^^. The idea is straightforward. 
We know that noninteracting composite fermions do not 
show FQHE at 5/2; our approach is to include the resid- 
ual interactions between them by constructing a basis 
of "noninteracting," or the "unperturbed," CF ground 
and excited states and rediagonalizing the Coulomb in- 
teraction in that subspace to obtain the spectrum for 
"interacting" composite fermions. This is known as the 
CF diagonalization, and the relevant techniques are de- 
scribed in the literature^^i^. As usual, we simulate the 
second LL physics in the lowest LL by working with an 
appropriate effective interaction. We work at the same 
flux value as the Pfaffian wave function, but because of a 
technical reason^^, we work with holes, rather than elec- 
trons. (Holes are not to be confused with quasiholes.) By 
particle hole symmetry, the number of holes is given by 
iV,, = {2Q + l)-N = iV-2 at 2Q = 2N -i. In what fol- 
lows, composite fermions are made by attaching vortices 
to holes rather than electrons. We show results at "ze- 
roth order" CF diagonalization (when only the lowest en- 
ergy unperturbed states are considered) and "first order" 
CF diagonalization (which also includes states with one 
higher unit of "kinetic energy"). The composite fermion 
kinetic energy levels are called A levels. 

Figure [To] shows the excitation spectra at the half filled 
second LL for Nh = 12, 14, 16 and 20 obtained by 
CF diagonalization at the zeroth and the first orders. 
[Nh — 18 is not considered as it aliases with v — ^ 
of holes.) The residual interaction between composite 
fermions lifts the degeneracy between various states to 
produce an incompressible state already at the lowest 
(zeroth) order, which neglects A-level mixing. Although 



the energy gaps change by up to 50% in going from the 
the zeroth to the first order, the incompressibility is pre- 
served, indicating that while A-level mixing renormalizes 
composite fermions, it does not cause any phase tran- 
sition. The overestimation of gaps at the zeroth order 
may be attributed to the very small dimensions of the 
CF basis. All CF basis states are perturbations of the 
noninteracting CF Fermi sea, making it explicit that a 
rearrangement of composite fermions near the CF Fermi 
level is responsible for the | FQHE. Although there is 
some ambiguity as to which excitation is to be identified 
with the transport gap (corresponding to a far separated 
quasiparticle-quasihole pair), an inspection indicates a 
gap of ~ 0.02, which is consistent with the earlier results 
from exact diagonalization^-'^^. 

Figure [TT] shows analogous results for the half filled 
lowest LL. The zeroth order CF diagonalization gener- 
ates the lowest band, and the first order generates the 
next band. The energies of states in the lowest band 
do not change appreciably from zeroth to the first order. 
The energy gap between the two lowest bands can be 
understood as the energy cost of exciting one more CF 
particle-hole pair. No such bands are seen for half filled 
second LL. 

It is not known how this description of the 5/2 FQHE 
relates to the Pfaffian model. In particular, a natural 
description of the quasiparticles is as excited composite 
fermions (which are heavily renormalized by the residual 
interaction). From this perspective, there is no reason 
to suspect that they would obey nonabelian statistics, 
although that cannot be ruled out as the residual inter- 
action causes nonperturbative change. 
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